A pulsed Sisyphus scheme for laser cooling of atomic (anti) hydrogen 
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We propose a laser cooling technique in which atoms are selectively excited to a dressed metastable 
state whose light shift and decay rate are spatially correlated for Sisyphus cooling. The case of 
cooling magnetically trapped (anti)hydrogen with the 1S-2S-3P transitions using pulsed ultra violet 
and continuous-wave visible lasers is numerically simulated. We find a number of appealing features 
including rapid 3-dimensional cooling from ~1 K to recoil-limited, millikelvin temperatures, as well 
as suppressed spin-flip loss and manageable photoionization loss. 

PACS numbers: 37.10.De, 67.63. Gh 



Recent progress [H HI in producing antihydrogen (H) 
improves the prospects for precision spectroscopy and 
de Broglie wave interferometry of H that may uncover 
new physics in low-energy experiments 0, Anti- 
hydrogen atoms are synthesized from anti-protons and 
positrons [IrlEI in such small numbers that trapping and 
efficient cooling to millikelvin temperatures are likely re- 
quired for precise measurements. In contrast to hydro- 
gen, which may be cooled by collisions with a buffer 
gas 0] , or by selecting low-energy atoms from an intense 
beam, cooling of H will likely rely on laser cooling tech- 
niques 

The most obvious approach to laser cooling of H Q 
or H demands 121.6 nm Ly-a radiation. Apart from the 
difficulty of manipulating VUV light, the generation of 
even ~10 nW of CW Ly-a radiation is technically chal- 
lenging Q. In addition, such cooling of H atoms in a 
kelvin-deep magnetic trap faces several interrelated diffi- 
culties. The need to avoid spin-flip losses, combined with 
the limited fraction of phase space addressable with a low 
intensity, single-frequency laser, implies that the cooling 
will be slow. Indeed, the only experimental Ly-a cooling 
work so far 0] used a pulsed laser with an average power 
of 160 nW (2.5 nW at the location of atoms) to cool mag- 
netically trapped H, and took more than 15 minutes to 
reach 8 mK, starting from just 80 mK. Furthermore, 3D 
cooling in this approach was aided by collisional mixing, 
which will be absent in dilute samples of H. Instead of re- 
lying on Ly-a radiation, there are several proposed cool- 
ing schemes using more readily available lasers to drive 
Doppler sensitive 2-photon transitions flol - |l~2j ] . However, 
in addition to limited phase-space addressability similar 
to Ly-a cooling, these schemes have the difficulty of losses 
due to photoionization. 

Motivated by previous work which used transi- 
tions between excited states for cooling and trapping, 
we propose a 3- level cooling scheme (Fig. [T|), applicable 
to magnetically trapped H, where a metastable state |e) 
is coupled to a short-lived state |e') by a blue detuned 
standing wave coupling f2 ee ' . Atoms in the ground state 
\g) are repeatedly excited to the bottom of the dissipa- 
tive |e) — |e') optical lattice by a pulsed, Doppler sensitive 
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FIG. 1: (a): Level diagram for the proposed cooling scheme 
(see text), (b): Simplified dressed-state picture of the cooling 
scheme. The pulsed two-photon excitation has bandwidth 
1/r, is detuned from the bottom of the lattice with depth U by 
So, and provides Doppler cooling by transferring momentum 
hk ge to an atom with velocity v. The linewidth T(x) of the 
dressed excited state \e(x)) is indicated by the width of the 
gray curve. As the atom climbs the hill the velocity (green 
arrow) decreases while the decay probability increases, leading 
to Sisyphus cooling. (The \e) —¥ \e) decay is ignored in (b). 
Not shown is the other, detuned dressed state \e'(x)).) 



2-photon coupling £l ge . The proposed cooling process 
arises from two effects: 2-photon Doppler cooling 
12j associated with \g) — > |e) excitation, and Sisyphus 
cooling [l4| associated with the |e) — |e') lattice. 

For magnetically trapped H, \g), |e) and |e') are the 
maximally Zeeman shifted states in the IS, 2S and 3P 
manifolds respectively. We show that the cooling scheme 
provides both a large capture velocity (100 m/s) and a 
low final temperature (near the Ly-/3 single-photon re- 
coil temperature of 1.8 mK), and allows for large vol- 
ume 3D cooling. Advantages include: availability of both 
nanosecond-pulsed UV 2-photon 1S-2S radiation [Iff and 
CW 2S-3P radiation at 656 nm; reduction of UV pho- 
toionization losses; and suppressed spin-flip transition 
from the 3P level in a high field. In the following we 
first discuss the pulsed Sisyphus cooling scheme in a ID, 
semiclassical model. After justifying the model with a 
ID quantum simulation fl6L Il7j. we present a 3D semi- 
classical simulation for magnetically trapped H. 

The proposed cooling scheme involves repeated pulsed 
excitations, each followed by spontaneous decay 18j. The 
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propagation direction of the 2-photon excitation pulses 
alternates between ±i. The Rabi frequency Q ge (x,t) is 
(0/r)f(t/T)e ±tk9 " x , where 9 is the pulse area, k ge is the 
sum of the wavevectors of the two photons, and f{t/r) 
is a normalized pulse-shape function with characteris- 
tic duration r. For 9 -C 1 its Fourier transform F(ujt) 
gives the excitation spectrum. The interval T rop between 
pulses is long enough so that excited atoms, moving in 
the |e) — |e') lattice, decay to \g) with high probability. 
In the effective 2-level system (Fig.[Tt>), the spatially de- 
pendent detuning 8{x) and lincwidth T(x) (see Eq. JT])) 
allows spatial selectivity in both pulsed excitation and 
subsequent decay. The |e) — \e') transition is driven by a 
standing wave coupling Q ee i (x) with a positive detuning 
A, resulting in two dressed states \e(x)} and \e'(x)) which 
are spatially dependent superpositions of |e) and |e') and 
connect to those states respectively as f2 ee < — > 0. The 2- 
photon \g) — |e) detuning from the unshifted mctastable 
state |e) (decay rate 7 es f=a 0) is 8q, and 7 e ' g ,7 e ' e are the 
decay rates from |e') (Fig. QJi). We assume A ^> S ,j e > g 
and 7 e 'g 3> 7 e ' e - Atoms are predominantly excited to, 
and adiabatically follow, the dressed state \e(x)) (and 
not \e'(x))) 0, and 
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In Eq. (TT]) we have ignored 7 e ' e and 7 eff , but their inclu- 
sion has little influence on the results described below. 
The depth U of the resulting |e) — \e') optical lattice is 
given by the maximum of Sq — S(x), and its period is 
determined by k ee i. We define kij as the wavevector of 
the i — j transition for i,j = g, e, e', with the associated 
recoil velocities and frequencies defined as Vij = Hkij/m 
and uj r sj — hkfj/(2m) where m is the atomic mass. 

Doppler-sensitive |o) — |e) absorption leads to 2-photon 
Dopplcr cooling [10Hl2j . The Doppler shifted, spatially 
dependent detuning, S(x, v) = 6(x) ± k ge v, allows atoms 
at different velocities to be excited to the lattice poten- 
tial at different locations. Atoms with velocity v are 
resonantly excited at positions such that |(5(x,w)|r < 1. 
For r > 2 7T / 1 <5 1 , strong excitation only occurs in the 
velocity range Vd < \v\ < v c , with decoupling velocity 
v d — \3o\/k ge and capture velocity v c ~ (\8q\ + U/h)/k ge . 
Atoms with \v\ <C Vd or \v\ ^> v c are off resonance and 
not efficiently excited. 

In addition to Doppler cooling, the correlation between 
the spatially dependent detuning 8(x) and the decay rate 
T(x) leads to Sisyphus cooling since atoms preferentially 
decay from |e(x)} at the tops of the light shift poten- 
tial I4J. The Sisyphus effect is particularly efficient for 
atoms with v close to Vd, which are excited near the bot- 
tom of the lattice. If |m» 2 < U, atoms remain within 
one lattice site and typically oscillate before decaying to 
\g). The decay is enhanced at the classical turning point, 
due to both the larger decay rate and the longer time 
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FIG. 2: OBE simulation of cooling properties. (a,d): Normal- 
ized excitation probability p e vs xo = k ee i xo/n and v = v/v ge . 
(b,e): Normalized energy loss per pulse e vs v. e is aver- 
aged over xo and is in units of Tiu) rige . (c,f): Ratio of e 
to the normalized ionization probability, s/pi vs v. Here 

5 = -25oj r , 9e , r = -2.5/<5 , /(*) = -^e-' 2 , 9 = tt/8, 
Q ee /(x) = f2 ee / sin(fc e e ,a 0> A = Q ee i /2 = 200 u) r , g e, and 
le'g = 87 e ' e = 2u) rige . The red curves in (b,c,e,f) correspond 
to Q, ce i = 0. The left and right panels are for k ge /k ee i =25 
and 5.4 respectively. 



spent there. Averaged over the position (and velocity) 
dependent decay probability, the velocity distribution af- 
ter decay is centered at zero velocity, with a rms width 
less than for v 3> v ge . On average this removes more 
than 75% of the atomic kinetic energy per 2-photon ex- 
citation. 

To characterize the cooling, we define the normal- 
ized position and velocity dependent excitation proba- 
bilities p e = 4P e /8 2 and energy loss per 2-photon pulse 
e = 4AE/6 2 , where P e is the 2-photon excitation prob- 
ability, AE is the energy loss per 2-photon pulse and p e 
and e are ^-independent for 9 <C 1. Figure 2(a,b,d,e) 
shows p e and e, determined using 3-level optical Bloch 
equations (OBE) for a "dragged atom" following a trajec- 
tory x(t) = x n +vt. Figure^a, b) and Fig.^d, e) are for 
k ge /kee' = 25 and 5.4 respectively. The latter ratio corre- 
sponds to fcis-2s/fc2S-3P in hydrogen. Figure [2] suggests 
that the simple picture of phase-space selective excitation 
plus Sisyphus cooling only applies for k ee ' <C k ge , (vt) , 
as in Fig. [2^a,b). For atoms that move more than l/k ee ' 
during r, the excitation is complicated by multi-photon 
resonances at velocities with k ge v + 2nk ee 'V rj 6 (the 
peaks of the black curve in Fig. [2^ for v > 20) . In 
addition, for large k ee ' , Doppleron resonant coupling to 
\e'(x)) [20j occurs at moderate speeds v with 2nk ee 'V ps A 
for integer n (the sharp dips of the black curve in Fig. [2^ 
near v = 70, 80), and leads to heating. Neverthe- 
less, efficient Sisyphus cooling is still possible for mod- 
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erate k ge /k ee > (see Fig. [2^). Compared to regular 2- 
photon cooling (red curves in Fig. [2jb,e)), the peak ex- 
citation probability is decreased by approximately ^jr/jr 
due to the spatially inhomogeneous broadening of |e(x)) 
(Fig. [I}}). However, due to the Sisyphus enhanced energy 
removal per excitation, the average energy loss per pulse 
remains comparable to the Doppler-only case, but with 
an increased velocity capture range Vd < \v\ < v c . 

In addition to the increased velocity capture range, 
the decreased excitation probability to |e(x)) also helps 
mitigate the photonionization loss from |e(x)) to the con- 
tinuum. For degenerate 2-photon excitation of H to the 
2S level, the ionization probability per pulse is given by 
/ dt^ ioni (t)p 2 s(t) where p 2 s{t) is the 2S state pop- 
ulation and 7i on i is the rate of ionization from 2S due to 
the UV radiation [2l[ . As a measure of cooling efficiency 
per 2-photon pulse, in Fig.|2fc,f) we compare e/pi, the ra- 
tio between the normalized energy loss e and normalized 
ionization probability pi — l2P lon i/0 3 , with (black curve) 
and without (red curve) the Sisyphus cooling. Here we 
have set 7i on j = 1.6f2 ge [2l[. We see that the Sisyphus 
effect enhances the cooling efficiency by approximately 
U/\h5o\ near v = Vd where the Doppler cooling has the 
best s/pi. 

We simulate the cooling process with a semiclassical 
stochastic wavefunction (SCSW) method lj| 17 1. The 
simulation of a cooling cycle is divided into two stages: 
excitation (0 < t < t) and decay (r < t < T rcp ) (we 
ignore quantum jumps during excitation). The exter- 
nal motion of the atom is described by a classical tra- 
jectory x(i). The internal dynamics arc described by 
a stochastic wavefunction \ip(t)), which, after the g — e 
pulse, is probabilistically projected to either \g) or the 
{|e), |e')} manifold (typically almost all in |e(x))) as 
\tp p ). Due to this post-selection, the optical force in 
the excitation stage cannot be evaluated in the usual 
way as (ip(t)\F\tf)(t)) where F is the force operator. In- 
stead, the force is estimated as the real part of a "weak 
value" [H, (,p p {t)\F\iP(t))/(iP P {t)\i;(t)), where |V(*)> is 
found by forward-propagating the pre-determined state 
IV'(O)) = \g) and is found by back-propagating 

the post-determined state (iP p (t)\, both for a dragged 
atom. This estimation method reproduces the quantum- 
mechanically expected velocity change during the pulse, 
impulse, due to both the recoil effect and the excited-state 
dipole force. During the second stage, the stochastic 
wavefunction \ip p (t)) € {|e),|e')} manifold and x(t) are 
propagated in small time-steps, until a quantum jump 
occurs 3, If the quantum jump is an |e') — > |e) tran- 
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FIG. 3: Comparison of ID SCSW and QSW simulations. Here 
n ee //(47r) = A/(2tt)=2.7 GHz as in Fig. H 6 = -2-k/t, 
9 — 7r/4. (a): Average speed (|«|) vs pulse number for three 
different pulse durations r. 4 ns represents a compromise 
between 1 ns, where the bandwidth is so large that there is 
little spatial selectivity, and 7 ns, where the bandwidth is so 
small that only a small fraction of atoms are excited. Thick 
lines are an average of 30 QSW trajectories, while thin lines 
are an average of 20 SCSW trajectories, (b): The equilibrium 
temperature T vs pulse duration r. At small velocities the 
SCSW method becomes less accurate for r > 9 ns, roughly 
set by half the oscillation period in a single |e) — |e'} lattice 
site. 



nal and external degrees of freedom of a 3-level atom, 
to confirm that the SCSW method correctly predicts 
the cooling dynamics and the final temperature. In 



Fig. E] typical results for smoothed square pulses 24 1 
are compared, for the appropriate hydrogen 1S-2S- 
3P parameters, r y e ' g /(2n) = 26.6 MHz, j e ' e /(2n) = 



3.6 MHz, {v ge , Wee' , Vge' } = {3.3,0.6,3.9} m/s, and 
{w r , ge ,w r , ee >,w r)ge >}/(27r) = {13.4,0.46,18.8} MHz, 
which are also used in Fig. HJd-f). We find good agree- 
ment between the SCSW and QSW methods as long as 
the dragged atom picture is valid during the pulse, i.e., if 
the optical force during the short excitation does not sig- 
nificantly displace the trajectory compared to the wave- 
length (k ee r Sv pu \ sc t ^ l)(Fig. Ep). The attainable ID 
temperature predicted by the quantum simulation de- 
creases with decreasing bandwidth 1/r, and is remark- 
ably low (~3 mK) even with t = 5 ns. 

Having verified the semiclassical approach for our pa- 
rameters, we use SCSW to simulate 3D cooling of mag- 
netically trapped H. To describe the 3D light-atom in- 
teraction, we included all ten electronic levels in the 
1S-2S-3P manifold (ignoring hyperfine structure). In a 
high magnetic field, the cooling process is dominated by 
the three maximally Zccman-shiftcd states of the IS, 
2S and 3P levels (corresponding to \g), |e) and |e')), 
which would form a closed system under 1S-2S 2-photon 
coupling Q ge and perfect a + 2S-3P coupling Q ee > . We 
consider a magnetic trap with B = {B x ,B y ,B z } = 



sition, we project \ip P (tY)to the dressed states \e{x)) or {Biy — B 2 zx/2, B\x — Bizyj2, B + B 2 (2z 2 



|e'(a;)) probabilistically [23J, while for an |e') — > \g) jump, 
we propagate x(t) freely until the next pulse. Upon each 
spontaneous emission, we use random velocity jumps to 
account for the recoil effect. 

We use a ID fu ll q uantum stochastic wavefunction 
(QSW) simulation |l6l. which includes both inter- 
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B Q = 0.75 T, Bi=0.8 T/cm and B 2 = 12 mT/cm 2 , sim- 
ilar to those for an existing antihydrogen apparatus 0. 
Both \g) and |e) feel a trapping potential V ~ Mb|B| 
(p-Q is the Bohr magneton), so the Uj) — |e) detuning 
8q is nearly free from Zeeman shifts [25J. The |e) — |e') 
detuning A(£?) w A — pbB/H, on the other hand, is 
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FIG. 4: Evolution of atomic velocity (a) and position (b) for 
a typical classical trajectory during the simulated cooling of 
magnetically trapped H. The insets give the quasi-equilibrium 
distributions. 



field-sensitive and has a position-dependent shift. 

We consider a 2S-3P lattice composed of three pairs of 
standing wave Gaussian beams, each with 1/e 2 diameter 
d, arranged symmetrically with equal intersection angles 
a to z. The choice of relative phases between standing 
waves is not critical to the cooling scheme. The beams 
are circularly polarized to maximize the rr + components 
(relative to B). In the Paschen-Back regime considered 
here, with fi B B > TiA 3Pfinc (A 3Pjfinc /(27r) = 3.25 GHz), 
the 7r coupling induces spin flip losses after 2S excitation 



with a branching ratio of r s f = 



/(A + ^B/ny 



3P,finc 

(similar for a~ coupling). Even for 2S-3P light that is 
purely 7r or a~ polarized, the spin-flip probability per 2S 
excitation is still less than 0.3% in a field of 1 T. 

Figure Q] plots a typical classical trajectory of H dur- 
ing cooling. The simulation starts with H in \g) at the 
trap bottom (B = 0.75 T), with initial longitudinal (z) 
and transverse (x,y) kinetic energy of Ei = 0.5 K and 
E t =0.25 K respectively. The 2-photon excitation beam 
overlaps with the 2S-3P beams in a 12 cm long and 1.8 
cm wide cooling zone, approximately covering the trap 
up to the 0.2 K equipotential surface (26[. We choose 
d = 3 cm and a = 0.1 for the 2S-3P beams, with peak 
intensity of 0.46 kW/cm 2 per beam corresponding to 
n ee ,/(2n) =1.3 GHz. A(0.75T)/(2tt) = 5.3 GHz is cho- 
sen so that A(B) ^> 7 e ' g within the cooling zone. In 
a flatter octopole trap [l[ the detuning constraint is re- 
duced, allowing for a reduced f2 ee < and less 656 nm power. 
The total power requirements can also be lessened with 
a moderate finesse optical cavity. We choose 9 = ir/8, 
t = 4 ns and T rep = 2 /xs [l^. The 2-photon detuning 
5q = —tt/t is chosen to improve the scattering rate for 
longitudinally cold atoms that are transversely hot. 

The rapid cooling trajectory shown in Fig. |4] is typi- 
cal for atoms with Ei < 0.5 K and E t < 0.25 K, which 
are cooled to quasi-equilibrium within A^otai = 4 x 10 4 
pulses during a cooling time of only 80 ms. While some 
atoms with E t significantly larger than 0.2 K may orbit 
around the cooling zone and not be efficiently cooled, the 
final velocity distribution for most atoms is remarkably 
isotropic with 5 m/s width (Fig. 0^ inset), which can be 
further reduced by increasing r (Fig.[3p). The total spin- 



flip loss is found to be less than 0.1%. As with other 
hydrogen cooling proposals fiol- 12|. one must consider 
limitations imposed by photoionization losses. We pcr- 
turbatively calculate photoionization from state popula- 
tions determined by internal state dynamics that ignore 
photoionization. For a reasonable two-color 2-photon ex- 
citation scheme where the stronger laser beam cannot 
ionize H from the 2S state in a single step [HI , we found 
ionization losses less than 25% |28| . Even if 243 nm radia- 



tion is used for the 1S-2S excitation [21| . photoionization 
loss would still be less than 25% for cooling, by applying 
iVtotai = 10 9 , 9 = 2.5 mrad pulses [H in 2000 s. 

We have proposed and analyzed a pulsed Sisyphus laser 
cooling scheme applicable to magnetically trapped H or 
H. The approach leads to rapid 3D cooling to <10 mK 
in a magnetic trap over a large volume with small spin- 
flip losses. Approaches to reduce photoionization losses 
to practical levels are proposed. Cooling efficiency may 
be further improved by exploring the spatial-temporal 
control of both the 2-photon excitation and the excited 
state lattice. This excited-state Sisyphus method may 
open new possibilities for cooling of deuterium, tritium 
or other species with metastable states. 

We would like to thank Amy Cassidy, Gretchen Camp- 
bell, and Jonathan Wrubel for helpful discussions. 



[9 
[10 

[11 
[12 

[13 
[14 
[15 



[16 
[17 

[18 
[19 



[20 
[21 
[22 
[23 
[24 



G. 
Y. 
M. 
G. 
G. 
R. 



B. Andresen et al, Nature 468, 673 (2010). 
Enomoto et al, Phys. Rev. Lett. 105, 243401 (2010). 
Amoretti et al, Nature 419, 456 (2002). 
Gabrielse et al, Phys. Rev. Lett. 89, 213401 (2002). 
Gabrielse et al, Phys. Rev. Lett. 100, 113001 (2008). 
deCarvalho et al, Can. J. Phys. 83, 293 (2005). 
Laser Cooling and Trapping, H. Metcalf and P. van 
der Straten (Springer- Verlag, 1999). 
I. Setija et al, Phys. Rev. Lett. 70, 2257 (1993). 
M. Scheid et al, Optics Express 17, 11274 (2009). 
M. Allegrini et al, Phys. Lett. A 172, 271 (1993). 
V. Zehnle et al, Phys. Rev. A 63, 021402 (2001). 
D. Kielpinski, Phys. Rev. A 73, 063407 (2006). 
S. Wu et al, Phys. Rev. Lett. 103, 173003 (2009). 
A. Aspect et al, Phys. Rev. Lett. 57, 1688 (1986). 
A two-color excitation scheme is proposed by L. Yatsenko 
et alia Phys. Rev. A 60, 4237 (1999). Also see T. Kanai et 
al, Opt. Exp. 17, 8696 (2009) for wavelength conversions. 
J. Dalibard et al, Phys. Rev. Lett. 68, 580 (1992). 
H. J. Carmichael, Lecture Notes in Physics, Vol. ml8 
(Springer, Berlin, 1993). 

T. Breeden et al, Phys. Rev. Lett. 47, 1726 (1981). 
For atom moving with velocity v, the adiabatic following 
requires (e(x)\vd x \e(x)) ~ k ee iv—j^- <C A. The excita- 
tion predominantly to \e(x)) requires k ge v <C A. 
J. J. Tollett et al, Phys. Rev. Lett. 65, 559 (1990). 
M. Haas et al, Phys. Rev. A 73, 052501 (2006). 
Y. Aharonov et al, Phys. Rev. Lett. 60, 1351 (1988). 
J. Chen et al, Phys. Rev. A 47, 2128 (1993). 
Rising/falling edges are smoothed by sine functions over 
20% of the pulse duration. 



5 



[25] C. L. Cesar et al, Phys. Rev. Lett. 77, 255 (1996). 

[26] To avoid excitation outside the cooling zone, one can 
overlap two excitation beams with distinct frequencies 
of cji and 012 = CJ2S — ui respectively, where L02S is the 
frequency of the 1S-2S transition. Also see ref. [TH |. 

[27] To allow slowest atoms confined by the |e) — |e') lattice 
to decay, the |e) — > |<?) decay rate (7 e9 ) is increased from 
to 10 /is" 1 for t > 1 /lis. This can be effectively realized 
with a broadband |e) — |e') coupling pulse. 



[28] The 2S-3P radiation also photoionizes H from 3P, which 
in the simulation leads to less than 4% total loss for cool- 
ing from Ei + Et at 0.75 K. Also see P. Antoine et al, J. 
Phys. B 29, 5367 (1996). 

[29] Reducing the pulse area while increasing the total pulse 
number suppresses the total photoionization loss, as the 
ionization loss oc 9 3 for 5 < 1. 



